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With the quantum interference between two transition pathways, we demonstrate a novel scheme 
to coherently control the momentum entanglement between a single atom and a single photon. The 
unavoidable disentanglement is also studied from the first principle, which indicates that the stably 
entangled atom-photon system with superhigh degree of entanglement may be realized with this 
scheme under certain conditions. 
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I. INTRODUCTION 

In recent years, entanglement with continuous vari- 
ables attracts substantial attention for its importance in 
quantum nonlocality [l| and quantum information pro- 
cessing (QIP) @. As a physical realization, momentum 
entanglement has been extensively studied both theoret- 
ically ]!, 0, 0, H, 0] and experimentally [1]. With momen- 
tum entanglement between atom and photon, it is possi- 
ble to define the best localized single-photon wavepacket 
even in free space and realize the highest degree of 
continuous entanglement Q up to date. 

As known, photon emitted from atom will recoil and 
be entangled with the atom [1, H, Q due to momentum 
conservation. In order to coherently manipulate the en- 
tanglement, in this paper, we propose a novel scheme 
to control the entanglement through the atomic spon- 
taneously generated coherence (SGC). With the config- 
uration in Fig. 1(a), we find that, the recoil-induced- 
cntanglement will be affected by the interference between 
different transitions from the atomic upper levels, and 
can be effectively controlled by an auxiliary coupling field 
if the dipoles for the transitions are parallel. This scheme, 
compared with the Raman scattering [3| and resonant 
scattering |6j , could be more efficient in producing super- 
high degree of entanglement, since the controlling light 
is classical and need not be far detuned. 

To consider the scheme in a more realistic situation, 
we study the unavoidable process of disentanglement fol- 
lowing the generation of entanglement. From the first 
principle, we obtain the master equation and the char- 
acteristic time scale for the disentanglement. Comparing 
the two processes, we yield an upper bound for the degree 
of entanglement that may be steadily produced with the 
scheme. Under realistic conditions [9(, it is shown that 
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FIG. 1: (a) Atomic configuration with two well-separated 
upper levels. The momentum entanglement can be controlled 
by the classical coupling field through the SGC. (b) Schematic 
diagram for the momentum detections. The detectors for the 
atom and photon are restricted in one-dimension which is 
perpendicular to the propagation direction of the coupling 
light. 

the robust atom-photon entangled pair can be produced 
with superhigh degree of entanglement. 

II. ENTANGLEMENT GENERATION BY 
INTERFERENCE 

Concerning the kinetic degrees of freedom, the Hamil- 
tonian for the system depicted in Fig. 1 (a) can be writ- 
ten with the rotating wave approximation (RWA) as: 

k i =1 
+hJ2 [gi{k)a cl a\e- fk - ? + g 2 {k)u c2 aie- i% - ? + H.C. 

k 

+h (n ie ~ Wot e l ^a lb + n 2 e- Wot e^a 2b + H.C.) , 
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where hp and r denote the atomic center-of-mass mo- 
mentum and position operators. b a p denotes the atomic 
operator |a)(/?| (a, (3 = 1,2,6, c), and (at) is the an- 



nihilation (creation) operator for the fcth photonic mode 
with wave vector fc and frequency luj: = ck, where we use 
k to include both the momentum and polarization of the 
photonic mode for simplicity. gj(k) is the coupling coeffi- 
cient for the transition \j) — > \c) (j = 1, 2) and ilj denotes 
the Rabi frequency for the coupling \j) «-> \b) (j = 1,2). 
For the convenience in further calculations, we may ne- 
glect the dependence on k for gj(k) and treat them as 
constants. Uq and fc denotes the frequency and the wave 
vector of the coupling light, and uj a p is used to represent 
the frequency difference as: uj a p = Lu a — uip. As the evo- 
lution is considered in a close system, it is convenient to 
expand the photon-atom state in the Schrodingcr picture 
as: 



[ai(q)W, 0, 1) + a 2 (q)\q, 0, 2) + b(q)\q, 0, b}} 
c(q,k)\q,k,c), (2) 



E 

q, k 



where the arguments in the kets denote, respectively, the 
wave vector of the atom, the photon, and the atomic 
internal state. 

From the Schrodinger equation one may yield the dy- 
namic equations for the system with the Born-Markov 
approximation. With the transformation to the slowly 
varying parts: 

a x {q) = e -™+«i« ^ ■ A x {q), (3) 

a 2 (q) = e -™+«^]* . A 2 (q), (4) 

b(q-k ) = e -m?-Zo)+v bc ]t . B ^ _ ^ (5) 

c(q,k) = e- l[T ^+^ ]t -C(q,k), (6) 



where T(p) = hp 1 /2m and Aj 
yield: 
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dt 
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1 dt 

■ dC(q,k) 
1 dt 



,— ^ = <L,<^B ( q-l, n }-^A 2 iq) 



Ai(q)e~ 
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- lA ^A 2 (q), (9) 



= g ie m<r>-T(q+k)+^ k -^]t . Ai ^ + q 
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where 71,2 denote the linewidthes for the two upper levels 
|1) and |2); and e = jl x ■ j2 2 /\jli\ ■ \jl 2 \ with j2j being the 
dipole moment for the transition \j) — > \c) (j — 1,2). 
As in the experiments [8[, we restrict the detections for 




FIG. 2: (a) Distribution of \C(Aq, Ak)\ 2 with 7 = 1, A = 
10, fi = 1, S = 0.01, 77=0.001. One sees that the central 
peak dominates the whole function, (b) Plot of the local 
amplification of (a), which shows that Aq and Afc are highly 
correlated in the central peak, (c) Plot of the function Q(Q) 
with A specified as 5, 15, 40, for the bold, thin and dashed 
fines, respectively, (d) Plot of O(A) with Q specified as 0.5, 
3, 10, for the bold, thin and dashed lines, respectively. 7 = 1. 



the photon and atom in one dimension, which is also 
perpendicular to the propagation of the coupling field, 
as depicted in Fig. 1 (b). 

In order to have strong interference between the tran- 
sitions |1) — > |c) and 1 2) — > |c), the dipoles for the tran- 
sitions should be parallel or antiparallel, i.e., e = ±1 
[13, EH; furthermore, the dressed states produced by the 
coupling field should be nearly degenerate, which may 
be fulfilled when Ai/A 2 ~ —71/72 [lOj. To be consistent 
with these restrictions, as in some experiments Q, we as- 
sume e = 1, 71 = 72 = 7, fii = O2 = fi; and A 2 = —A, 
Ai = (1+5) A, where S is a dimensionless small term con- 
trolled by the detuning of the coupling field. The atom 
is initially prepared in state \b) with momentum wave- 
function as G(q) oc exp [— (q/5 p ) 2 ], where S p denotes its 
momentum variance. 

With the above simplifications, from Eqs. (7) to (10), 
it is straightforward to yield the solutions for the whole 
system. As a result, the steady state atom-photon en- 
tangled wave function reads: 



C(q,k,t-> 00) « xo- 



exp [-(Aq/r]) 2 } 
-Ai/7 + i(Aq + Afc)' 



(11) 



The Aq, Afc rj and Ai are defined in the appendix, where 
the other related mathematical details are also given, xo 
is a normalization factor. 

Theoretically, the entanglement of a pure state bipar- 
tite system can be completely evaluated by the Schmidt 
number K [12], which is defined as an estimation of the 
number of modes that make up the Schmidt decomposi- 
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FIG. 3: Relations between the Schmidt number K and the 
coupling field with 7=1, u X2 = 20, 77 = 0.001. (a) The bold, 
thin and dashed lines are plotted with Q = 0.1, £7 = 2 and 
f2 — 4, respectively, (b) The bold, thin and dashed lines are 
plotted with A = 9.97, A = 9.98 and A = 9.99, respectively. 



degree of entanglement K with respect to the detuning 
A = ^0 — <^2b and Rabi frequency Q, of the coupling field. 
From Fig. 3 and Eq. (12), one sees that the entanglement 
can be effectively enhanced by controlling the coupling 
field, and when v$ — ► LU2b + ^12/2, we have K — ► 00. 

Different to some previous works [1, [f|, the super- 
high entanglement produced with this scheme is due to 
the atomic spontaneously generated coherence (SGC) 
[ToL [ill ]. With parallel dipoles (e = 1), the photon 
emitted along the two transition pathes |1) — > |c) and 
|2) — > |c) will interfere and modify the momentum en- 
tanglement with the recoiled atom as a result. In recent 
studies, it is found that, with nearly degenerate upper 
levels and proper atomic coherence, the atom with SGC 
may exhibit anomalous enhancement of momentum en- 
tanglement in the spontaneous emission process (Tlj . In 
this scheme, however, it is shown that, even with well- 
separated upper levels and realistic conditions [§| , the en- 
tanglement could also be highly increased with the quan- 
tum interference which can be controlled by a classical 
light field. Therefore, the proposed scheme can most 
probably be used to produce highly entangled atom- 
photon pairs in realistic applications. 



tion. From Eq. (11), it is found that 0,0] 



K 



1 + 0.28 



|Re(A 1 / 7 )| 
0.28t? 



|9(A/ 7 ,rV7)l<5 2 ' 



(12) 



where the function 0(A/7, is defined in the ap- 

pendix and depicted in Figs. 2 (c) and (dV With respect 
to the realistic detections in experiment [8J], however, the 
degree of entanglement is better to be characterized by 
the ratio (i?) of the unconditional and conditional vari- 
ances in the momentum detections [ESQ, i.e.: 



_com 



R = 5q sinslc /Sq 
jingle = ^2 

Sq coia = (f) coin - (J>) 2 coln , 



(13) 
(14) 

(15) 



where <5q smglc is defined as the variance of the atomic 
momentum with single-particle detection, and 5q cmn de- 
notes the variance obtained by the coincidence detection 
on both the atom and the photon. From Eqs. (11) and 
(13), we yield: 



R 



V 



1.6|Re(Ai/ 7 )| 1.6|e(A/ 7 ,n/ 7 )|<S 2 



2.2K,(16) 



where one sees that the Schmidt number K can therefore 
be obtained experimentally by measuring the R ratio [1] . 

With Eq. (12) or (16), one sees that the entanglement 
is most sensitive with the relative detuning 5, therefore 
with the frequency of the coupling field. With the rela- 
tion UJ12 — Ai — A2 = (2 + 6) A, we plot in Fig. 3 the 



III. DISENTANGLEMENT 



As in the preceding section, to study the generation of 
momentum entanglement [III [Si, it is usually conve- 
nient to assume the entangled system to be a close pure 
state system. However, in a realistic environment with 
T 3> OK, the interaction with environment will make the 
entangled system into a mixed state, and, as a result, 
cause the disentanglement. Actually, only when the dis- 
entangling process is much slower than the generation 
of entanglement, the entangled system can well be ap- 
proximated by a pure state. With these considerations, 
it is then possible to give out an upper bound for the 
entanglement that could be produced reliably in the en- 
vironment. 

Concerning the momentum entanglement, the disen- 
tanglement is caused by the momenta exchange with 
the environment which may be composed of background 
atoms and photons. Theoretically, the influence from 
the background atoms can be eliminated by using a high 
vacuum system; therefore, in order to study the unavoid- 
able disentanglement, we can simplify the environment 
as a heatbath of background photons, coupled only to 
the entangled atom, as shown in Fig. 4 (a). In order 
to give a general analysis for this incoherent process, the 
atom is simplified as a two-level system with resonant 
frequency uj a , then the Hamiltonian of the total system 
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under RWA is: 
H to t = 
H s -- 

Hi = 



Hg + Hb 

(hp) 2 

2m 



Hi 



hw a &22 + huj k b\b k 



k 

hY, [g(k)a 12 ale- ikr 



H.C. 



(17) 
(18) 

(19) 

(20) 



where Hs, Hb and H\ denote the Hamiltonians for the 
system (the entangled atom-photon pair), the heatbath 
and the interaction between them, respectively. b k (b\.) 
is the annihilation (creation) operator for the entangled 
single-photon in its fcth mode, whereas a k and a k are 



those for the photons in the heatbath. 

It is known that the density matrix of the total system 
Ptot obeys the Liouville equation: 



Ptot = Aot(Ptot) = (£s + C B + A)pt 



(21) 



where the superoperators £tot j £s > £b and L\ are defined 
as: £tot(*) = — \[H\nt->*]i etc In order to reveal the 
dynamic evolution for the entangled system, we should 
adiabatically eliminate the heatbath terms from Eq. (21) 
to obtain the master equation for the entangled system. 

To proceed, we define the reduced density matrix for 
the system as p = Tr B (ptot) and a "projection state" as 
v = Tre(ptot) <8> Pb, where the trace "Tre" is taken over 
the heatbath space and pb denotes the initial state of 
the heatbath. As the coupling is weak, from Eq. (21), 
we yield the equation for the projection state as: 



v = Csv - ipB <8> Tr B ( £i J dr 



^.g(fc)a[e-™ fcT A(r)+H.C. , p® p B 



L k 



r 



(22) 



where the Markov approximation and the nonrelativis- From Eq. (22), we yield the master equation with the 
tic approximation hk/mc -C 1 are used, and A(r) = Lindblad form [13l ] as: 



e —iHsT/Hfi. e -ikr e iH S T/H 



P = ~bi H S, p] - T a [D(a n p + pan) + (1 + D)(a 22 p + po~ 22 ) 

D"P \ * „ik a r-~ „ — ik a r , ~ „—ik a r ~ ik a r~\ 

+ l>Y a |02ie po\ 2 e. +cr 2 ie p<J\ 2 e. J, 



(1 + D)T a [a 12 e- lKr pa 21 e^ r + a 12 e* k " r pa 21 t 



-ik a r~\ 



(23) 



where D = TrB(fit a ka ps) is the average number of the 
resonant photons in the heatbath; the atomic linewidth 
is given as T a = J2 k 2ir\g(k)\ 2 5(uj — oj a ). It is natural to 
assume that the heatbath is initially in the thermal equi- 
librium, i.e., pb = e" ffB/fcBT /Tr B (e" HB/A;BT ), then we 
have D = \j{^o.lk B T _ w here T is the temperature 
of the heatbath. 

The master equation, Eq. (23), describes the process 
of disentanglement, where the entangled system is now 
in a mixed state due to the interaction with the heat- 
bath. Theoretically, the entanglement of a mixed bipar- 
tite system can be evaluated with the "entanglement of 
formation" 14] . However, in order to base our analysis on 
direct experimental test @, we use the defined R ratio 
[cf. Eq. (13)] as the evaluation of the entanglement. 

As in Eqs. (14) and (15), the momentum variance of 
the single-particle (the atom) measurement is calculated 



while that of the coincidence measurement is: 
<2 r- 2, 



Sq c 



^single = r^p) _ [Tr^p)]*, 

icidence measureme 
Ei=i/ d « P(l> k o,hq,k ,i) 



(24) 



(25) 



where p(q,k,i,q',k',i') denotes the matrix element 
(q,k,i\p\q',k',i'), and the photon is assumed to be de- 
tected with some momentum kg. 

With Eqs. (23) to (25), it is straightforward to get: 



° Sq^/ADkl + RoTat' 



(26) 
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equality: 



2x10 s 4xl0 6 6x10 s 8x10 s 1 x 10 7 
( I. 



FIG. 4: (a) Schematic diagram of the disentanglement. The 
atom is treated as a two-level system and the environment 
is simplified as a heatbath of photons which is coupled to 
the atom through momenta exchange, (b) R(t) is plotted 
with Sq s inglc = 16kl and _R =3000. The temperature of the 
environment is specified as T — 300 K for the bold line, T = 
270 K for the thin line and T = 250 K for the dashed line. 



where Ro and <5qQ lnglc are denned as their initial values, 
Ro = R(t = 0) and Sq & iDglc = Sq sinsie (t = 0). It can be 
seen from Fig. 4 (b) that, R(t) decreases monotonously 
with time, therefore, the characteristic time scale Atdis 
for the disentanglement can be defined as R(t = Atdis) = 
\Ro. As Rq 3> 1, we yield: 



At, 



8q s insie 
ADRoklTa • 



(27) 



From Eq. (27) , one sees that the "disentangling time" 
Atdis is inversely proportional to the average number of 
the resonant photons in the heatbath [D{u> a )]. There- 
fore, by decreasing the temperature of the environment 
it is possible to significantly increase Atdis and then 
make the entangled system quite robust in the environ- 
ment. However, as the temperature can never reach ab- 
solute zero, this kind of disentanglement is "unavoid- 
able". Furthermore, as in Eq. (27), the disentangling 
time is also dependent on the initial entanglement, i.e., 
Atdis oc 1/Ro, which indicates that, the better entan- 
gled system is more fragile in the environment. Since 
Atdis — > when Ro — > oo, the ideal continuous EPR 
state [l[ can never be reached in a realistic environment 
in this sense. 

The physical meaning for the dependence on k a and 
T a in Eq. (27) is apparent: with larger energy and 
shorter lifetime for the transitions, the environment will 
exchange more momenta with the entangled atom per 
unit time, and as a result, accelerate the disentangle- 
ment. When temperature is low, the denominator in Eq. 
(27) has a sharp peak for the coupled frequency uj a , which 
ensures the two-level approximation reasonable for our 
treatment. 

As stated at the beginning of this section, the entangle- 
ment generated with Eqs. (12) and (16) is applicable only 
if the disentangling time is much longer than the time 
scale for producing the entanglement, i.e., Atdis 3> At cnt . 
Therefore, with the Eqs. (27) and (A8), we yield the in- 



R < 0.2* 



/ hk lc 5l 
DmklY a 



(28) 



which gives an upper bound for the entanglement that 
can be produced in realistic environment with this 
scheme. 

As in some reported experiments the atomic con- 
figuration with SGC as in Fig. 1 (a) can be realized 
by sodium dimers. With the experimental conditions 
7 ~ 10~ 7 wi c , oj\2 ~ IO7, f2 ~ 7, when the coupling 
field is tuned to A ~ W12/2 ~ 57 with S — 10~ 2 , from 
Eq. (A8), we have At cnt ~ 1ms. Take the time-of-flight 
into account the initial momentum variance can be 
prepared as hS p /m = lm/s, and then from Eq. (16) we 
obtain a superhigh degree of entanglement as R w 4600 
and K « 2100. To consider the disentanglement, we 
take Lo a = 5 x 10 14 Hz and r a = 10 7 Hz for estimations. 
With the environment temperature T — 150K, from Eq. 
(27), we have D ~ 10 -11 and Atdis ~ 10 4 s. One sees 
that the relation At ont <§; Atdis can be well fulfilled. 
Therefore, under these conditions, the robust highly en- 
tangled atom-photon pairs can be steadily produced in 
the environment. Actually, from Eq. (28), we have an 
upper bound as R <C 10 7 , which implies a strong abil- 
ity of producing entanglement with this scheme. On the 
other hand, if the environment is at a high temperature, 
e.g., T = 400 K, the disentanglement will be strongly 
enhanced and we now have At cnt ~ Atdis- With direct 
detections [8j , it is then possible to observe all these phe- 
nomena in experiment. 



IV. CONCLUSION 

In this paper, we demonstrate a novel scheme to pro- 
duce superhigh momentum entanglement between a sin- 
gle atom and a single photon with the atomic SGC [111 ]. 
Under certain experimental conditions @ , we show that 
the entanglement can be effectively controlled by the clas- 
sical coupling field and may be very robust against the 
disentanglement due to the environment. As we analyze 
the two physical processes separately and both from the 
first principle, most of our conclusions can directly apply 
to the previous models [3, |j, [5|, ly] . 

To give a better upper bound than Eq. (28), one may 
consider the generation of entanglement together with 
the disentanglement at the same time, which is also 
necessary to analyze the system when At cnt ~ Atdis 
rigorously. However, this method is more complicated to 
be generalized and will not change our above conclusions 
qualitatively. We plan to give the details of this method 
elsewhere in a future work. 
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Appendix: steady solutions of equations (7) to (10) 

In order to obtain Eq. (11), we define a matrix M as: 



M = 



£V7l72 



c \/TiT2 



in;, 



in* 



(Al) 



its eigenvalues and eigenvectors are denoted, respectively, 
as \j and (a,- , j3j , £j) T with j = 1,2,3. Then from Eqs. 
(7) to (10), the steady state solution of the entangled 
wave function can be written as: 



C(q, k, t — > oo) 



(A2) 



E 

3=1 



i(giotj + g 2 Pj)pjG(q + k-k ) 



-\j + % 



T(q) - T(q + k) + uj^ - (uj bc + v ) 



where pj is determined by the initial conditions, and 
when the atom is initially in \b), the restrictions are: 

3 3 3 

5>^=0, = E^0 = 1- (A3) 

3 = 1 3 = 1 3 = 1 

As the detection is restricted in one dimension, the solu- 
tion can be simplified as: 



C iq ,k)o,±^ + ^ e - {Aq/ny 

3 = 1 
3 



-Xj/j 2 + i{Aq + Ak) 



3 = 1 



(A4) 
(A5) 



where the effective wave vectors are defined as: 



Aq 



_ h(uj bc 



V<3) , Ubc 

-(q 



Ak 



mcj2 
ck - (u bc + v ) 

72 



(A6) 
(A7) 



and r\ = h{uj bc + v^)5 v /m^2C- We use Li ) 2,3(Ag, Ak) 
here to denote the three different terms that make up 
the summation in Eq. (A4). From Eq. (A4), it can be 
proved that: 



dA 9 dA/c|L 2l 3(A<?, Ak)\ 2 = o(n 2 /7 2 ) when n 0, 



which indicates that the L\(Aq, Ak) dominates the sum- 
mation when the coupling field is weak, as shown in Figs. 
2 (a) and 2 (b); therefore, Eq. (A4) can be well approxi- 
mated by the single-peak function as in Eq. (11). 

To give further analysis for the entanglement, from 
Eq. (Al), we yield Re(Ai/7) ~ ©(A/7, n/ 7 )5 2 , where 
the value of 0(A/7, n/7) is of order 0.1 or smaller as 
shown in Figs. 2 (c) and (d). Moreover, the time scale 
for producing the entanglement can be characterized as: 
At cnt = l/|Re(Ai)| = 1/7|6|5 2 , with Eq. (16), it may be 
written as: 



AU 



1.6R 



(A8) 



[1] A. Einstein, B. Podolsky, N. Rosen, Phys. Rev. 47, 777 
(1935); J. C. Howell et. al., Phys. Rev. Lett. 92, 210403 
(2004). 

[2] S. L. Braunstein and P. V. Loock, Rev. Mod. Phys. 77, 
513 (2005). 

[3] K. W. Chan, C. K. Law, and J. H. Eberly, Phys. Rev. 

Lett. 88, 100402 (2002). 
[4] K. W. Chan et. al, Phys. Rev. A 68, 022110 (2003); J. 

H. Eberly, K. W. Chan and C. K. Law, Phil. Trans. R. 

Soc. Lond. A 361, 1519 (2003). 
[5] M. V. Fedorov et. al, Phys. Rev. A 72, 032110 (2005). 
[6] R. Guo and H. Guo, Phys. Rev. A 73, 012103 (2006). 
[7] M. V. Fedorov et al, Phys. Rev. A 69, 052117 (2004). 
[8] M. D. Reid and P. D. Drummond, Phys. Rev. Lett. 60, 

2731 (1988); Michael S. Chapman et al, Phys. Rev. Lett. 

75, 3783 (1995); Christian Kurtsiefer et al, Phys. Rev. 

A 55, R2539 (1997). 



[9] H. R. Xia, C. Y. Ye, and S. Y. Zhu, Phys. Rev. Lett. 77, 
1032 (1996). 

[10] S. Y. Zhu and M. O. Scully, Phys. Rev. Lett 76, 388 
(1996). 

[11] arXiv: R. Guo and H. Guo, |quant-ph/0701018'V 2; S. Y. 
Zhu, R. C. F. Chan and C. P. Lee, Phys. Rev. A 52, 710 
(1995). 

[12] R. Grobe et al, J. Phys. B 27, L503 (1994); S. Parker 
et al, Phys. Rev. A 61, 032305 (2000); C. K. Law, I. A. 
Walmsley, and J. H. Eberly, Phys. Rev. Lett. 84, 5304 
(2000); W. C. Liu, S. L. Haan, R. Grobe, J. H. Eberly, 
Phys. Rev. Lett. 83, 520 (1999). 

[13] G. Lindblad, Comm. Math. Phys. 48, 119 (1976). 

[14] C. H. Bennett, D. P. DiVincenzo, J. A. Smolin and W. 
K. Wootters, Phys. Rev. A 54, 3824 (1996). 



